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>  The  dissociative  recombination  of  OPCl+  and  OPCl2+  has  been  studied  at  the  storage  ring 

Q_  CRYRING.  The  rate  constants  as  a  function  of  electron  temperature  have  been  derived  to  be 

Q  7.63  X  10_7(Te./300)'°89  and  >  1.2  X  10~6(Ttf/ 300)-1  22  cm3  s”1,  respectively.  The  lower  limit  quoted 

O  for  the  latter  rate  constant  reflects  the  experimental  inability  to  detect  all  of  the  reaction  products. 

The  branching  fractions  from  the  reaction  have  been  measured  for  OPCl+  at  ~0  cV  interaction 
energy  and  are  determined  to  be  MO+P+Cl)  =  (  16  ±  7)%,  /V(O+PC0  =  (!6±  3)%  and  MOP 
^  +C1)  =  (68  ±  5)%,  These  values  have  been  obtained  assuming  that  the  rearrangement  channel 

forming  P  +  CIO  is  negligible,  and  ah  initio  calculations  using  GAUSSIAN03  are  presented  for  the  ion 
structures  and  energetics  to  support  such  an  assumption.  Finally,  the  limitations  to  using  heavy  ion 
storage  rings  such  as  CRYRING  for  studies  into  the  dissociative  recombination  of  large  singly 
charged  molecular  ions  are  discussed.  ©  2008  American  Institute  of  Physics. 
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I.  INTRODUCTION 

Much  recent  environmental  debate  has  focused  on  the 
role  that  motor  vehicles  play  in  everyday  life,  ranging  from 
the  use  of  biofuels  and  electric  cars  to  the  efficiency  of  fuel- 
air  combustion  engines.  Several  aspects  related  to  the  latter 
of  these  areas  have  been  to  look  at  ways  of  making  the  fuel 
itself  cleaner,  e.g.,  by  including  additives  which  can  either 
control  or  effect  the  amount  of  environmentally  sensitive 
chemicals  released  from  the  combustion  process — in  the 
same  way  that  a  catalytic  converter  cleans  up  the 
exhaust — or  to  increase  the  efficiency  of  the  combustion  pro¬ 
cess,  which  would  then  have  the  knock  on  effect  that  the 
energy  output  from  the  process  could  be  maintained  but  that 
less  fuel  is  consumed.  Ion-molecule  chemistry  has  been 
shown  to  increase  the  rate  of  hydrocarbon  combustion  when 
ions  are  present  in  the  combustion  environment1  4  while 
oxyphosphorus  compounds  have  been  shown  to  be  effective 
replacements  for  phosphine  fuel  additives.5  Consequently, 
several  studies  on  the  ion-molecule  chemistry  of  phosphorus 
oxychloride  compounds  have  been  done  using  flow  tube  in¬ 
struments  in  order  to  understand  the  transition  from  the  pri¬ 
mary  electron  attachment  products  of  OPCl3  to  the  stable, 
fully  oxygenated  anion  species,  P03“  and  P02~  (Refs.6-10). 
These  two  ions  have  been  observed  in  the  flame  environment 
to  be  the  primary  terminal  ions.1112 

One  of  the  possible  competing  reactions  for  removing 
electrons  from  the  plasma  is  dissociative  recombination 
(DR).  DR  involves  the  attachment  of  a  typically  low  energy 
(<1  eV)  electron  to  a  singly  charged  molecular  ion  which 
then  dissociates  into  neutral  fragments.  Indeed,  DR  has  also 


a)Eleclronic  mail:  vz@physto.se. 


been  suggested  as  a  possible  reaction  process  that  increases 
the  efficiency  of  combustion  engines  since  it  leads  to  the 
production  of  reactive  radical  species  which  can  then  react 
further.13  Several  techniques  have  traditionally  been  used  to 
measure  the  DR  reaction  rate,  as  well  as  the  products  that 
arise  from  the  reaction,  from  flow  tube  instruments  to 
merged  beams  experiments.  One  extension  of  the  latter  tech¬ 
nique,  which  has  proven  to  be  especially  successful  in  the 
study  of  DR,  has  been  the  heavy  ion  storage  ring  since  it 
combines  the  ability  to  measure  to  high  precision  both  the 
reaction  rate  constant  and  the  reaction  products  with  the 
same  experimental  setup,  something  which  is  not  usually 
possible  in  flow  tube  instruments.  Over  the  last  ten  years  a 
considerable  body  of  data,  i.e.,  cross  sections,  rate  constants, 
and  reaction  products  has  been  generated  from  DR  experi¬ 
ments  at  both  storage  rings  and  flow  tube  instruments,  and 
several  reviews  of  these  data  and  comparisons  between  data 
sets  have  recently  been  published.14  16 

In  spite  of  the  uniqueness  of  the  information,  the  current 
generation  of  ion  storage  rings  has  several  limitations,  in¬ 
cluding  the  ability  to  conduct  experiments  at  only  300  K 
and,  of  relevance  to  DR,  for  only  relatively  low  mass  singly 
charged  molecular  ions.  The  former  places  some  constraints 
on  the  asymptotic  internal  temperature  of  the  ions,  especially 
if  their  internal  modes  couple  effectively  with  the  300  K 
black-body  radiation.  The  latter  arises  from  the  requirement 
to  velocity  match  the  ion  and  electron  beams  for  ~0  eV 
collisions,  which  is  experimentally  determined  by  the  ability 
to  create  and  confine  a  sufficiently  dense  electron  beam  to 
allow  the  reaction  to  proceed.  In  CRYRING,  the  ion  energy 
in  the  laboratory  frame  is  given  by  96 q2! M,  where  q  and  M 
are  the  charge  state  and  mass  of  the  ion,  respectively.  For  a 
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singly  charged  molecular  ion  with  mass  of  96,  this  corre¬ 
sponds  to  1  MeV  total  beam  energy  and  gives  rise  to  an 
electron  beam  of  only  5.7  eV,  The  low  kinetic  energy  of  the 
ions  also  places  constraints  on  the  unique  identification  of 
the  reaction  products.  The  detectors  used  in  storage  rings  to 
determine  the  products  from  the  DR  reaction  are  ion- 
implanted  silicon  detectors  which  essentially  convert  the  ki¬ 
netic  energy  given  up  by  the  particles  in  the  active  volume  of 
the  detector — through  the  creation  of  electron-hole  pairs — to 
a  proportional  output  voltage.  The  resolution  is  determined 
by  the  electrical  noise  generated  in  the  crystal,  and  if  the 
noise  is  the  same  order  as  the  signals  arising  from  real  par¬ 
ticles,  then  unique  identification  of  the  reaction  products  be¬ 
comes  difficult.  Some  solutions  to  this  problem  have  been 
suggested,  but  their  range  of  applicability  is  limited,  and 
these  are  also  discussed  in  the  recent  review  articles.14  lb 

The  present  study  addresses  both  of  the  issues  raised 
above:  an  investigation  into  the  DR  of  heavy  singly  charged 
phosphorus  oxychloride  compounds  of  relevance  to  combus¬ 
tion  additives,  e.g.,  OPCl+  and  0PC12\  The  DR  of  0PC1+  at 
0  eV  interaction  energy  can  proceed  through  four  available 
reaction  channels 


OPCT  +  e 


P  +  O  +  CI+  1.4  eV 

(a). 

O  +  PCI  +  4.6  eV 

(b). 

OP  +  Cl +  7.5  eV 

(c), 

P  + CIO +  4.2  eV 

(d). 

(0 


It  is  noted  that  channel  (Id)  involves  the  recombination  of 
the  two  “outside”  atoms  as  determined  from  the  structure 
calculations  which  are  discussed  later  In  the  DR  of  OPCl2+, 
the  following  channels  are  energetically  available: 

rOPCl  +  Cl  +  7.2  eV  (a), 

OP  +  Cl2  +  5.9  eV  (b), 

OP  +  Cl  +  Cl +  3.4  eV  (c), 


OPCl2+  +  ?  -v  i 


PCI  +  CIO  +  3.3  eV 


P  +  C120+  1.4  eV 


(d) , 

(e) , 


(2) 


PCI  +  Cl  +  0  +  0.5  eV  (f), 


^PC12  +  O  +  3.85  eV  (g). 


The  exothermicity  of  these  reaction  channels  [Eqs.  (1)  and 
(2)]  has  been  determined  using  both  Ref.  17  and  the  calcu¬ 
lations  presented  later  in  this  paper.  Due  to  the  relatively 
high  mass  of  the  ions,  some  difficulties  have  been  encoun 
tered  putting  constraints  on  the  measurement  of  both  the  re¬ 
action  rate  constant  and  the  reaction  products,  i.e.,  it  has  not 
been  feasible  to  measure  the  complicated  branching  fractions 
in  the  DR  reaction  of  0PC12+. 


II.  EXPERIMENT  AND  AB  INITIO  CALCULATIONS 


The  measurements  were  made  at  the  CRYRING  facility 
located  at  Stockholm  University.  The  apparatus  has  been  de¬ 
scribed  in  detail  previously18  and  so  only  the  salient  details 
which  are  relevant  to  the  present  measurements  will  be  men¬ 
tioned  here.  Both  0PC1+  and  0PC12+  were  created  from 
OPCI3  in  a  hot  filament  ion  source.  In  order  to  save  the  liquid 
and  the  filament  from  the  corrosive  nature  of  OPCI3,  the 


vapor  was  pulsed  into  the  source  in  conjunction  with  the 
injection  cycle.  In  doing  so,  the  high  concentration  of  OPCI3 
necessary  for  the  production  of  a  usable  amount  of  ions  was 
only  present  in  the  source  for  short  periods  of  time.  The  ion 
of  choice  was  mass  selected,  injected  into  the  ring,  and  ac¬ 
celerated  to  the  maximum  allowed  energy.  The  isotopologue 
OP,7CT  was  used  since  it  gives  a  slightly  better  mass  reso¬ 
lution  when  determining  the  reaction  products.  Furthermore, 
OP37Cl+  is  isobaric  with  S4Kr+  and  since  it  is  easy  to  create 
large  ion  currents  with  the  latter  ion  this  simplified  the  set¬ 
ting  up  of  the  experimental  ring  parameters,  i.e.,  the  ion 
beam  transport  and  storage  in  the  ring,  and  made  switching 
to  OPCI+  relatively  easy  in  spite  of  the  relatively  low  flux. 
Since  product  determination  was  impossible  for  the  heavier 
species,  the  dominant  isotopologue  OP35Cl2+  was  used.  Once 
in  the  nng,  the  ions  were  subjected  to  the  normal  measure¬ 
ment  cycle  including  a  period  in  which  the  ion  beam  was 
cooled  in  the  electron  cooler.14  Electron  cooling  implies  a 
reduction  of  the  phase-space  occupied  by  the  ions  through 
their  interaction  with  the  cold  electron  beam  characterized  by 
longitudinal  and  transverse  energy  spreads  of  0.1  and  2  meV, 
respectively.  Under  the  cooling  conditions,  the  electron  and 
ion  beams  have  the  same  velocity,  i.e.,  have  a  relative  colli¬ 
sion  energy  close  to  0  eV.  However,  due  to  the  heavy  nature 
of  the  10ns,  the  cooling  period  was  not  long  enough  to  ensure 
that  the  ions  reached  the  same  energy  spread  as  the  electrons. 
As  such,  the  main  purpose  of  the  electron  cooler  in  the 
present  experiment  was  only  to  serve  as  an  electron  target. 
However,  both  ions  have  nonzero  dipole  moments  and  the 
cooling  period  was  sufficient  for  spontaneous  photon  emis¬ 
sion  leading  to  vibronically  ground  state  ions  The  electron 
cooler  cathode  voltage  was  tuned  away  from  the  cooling 
conditions  which  enabled  us  to  vary  the  relative  collision 
energy.  The  cathode  voltage  was  ramped  in  such  a  way  that 
the  electrons  were  initially  faster  than  the  ions  and  then  their 
relative  velocity  decreased  until  they  were  slower.  As  such, 
the  relative  collision  energy  varied  from  1  eV 
through  0  to  1  eV  again.  Monitoring  the  count  rate  of  the 
DR  products  dNIdt  as  a  function  of  the  cathode  voltage,  the 
energy-dependent  cross  section  was  calculated  as 


oi,E)  - 


1  dN  V'V^-nr2' 

dt  /,/,•/ 


(3) 


where  t/rd  is  the  relative  velocity  between  electrons  and  10ns, 
ve  and  Vj  are  the  electron  and  ion  velocities  in  the  laboratory 
frame,  respectively,  le  and  /,  are  the  electron  and  ion  cur¬ 
rents,  respectively,  e  is  the  elementary  charge,  re  is  the  elec¬ 
tron  beam  radius,  and  /  is  the  length  of  the  electron  cooler. 
The  obtained  result  was  corrected  for  the  toroidal  regions  of 
the  electron  cooler  and  space-charge  effects  induced  by  the 
electrons  2°  1 

The  heavy  nature  of  the  ions,  at  84  and  1 17  amu,  results 
in  the  relatively  low  kinetic  energies  of  the  beams,  1.12  and 
0.8  MeV,  respectively.  This  led  to  several  experimental  dif¬ 
ficulties,  c.g.,  due  to  the  insufficient  electron  cooling,  the  ion 
beam  profiles  were  quite  broad  which  could  lead  to  products 
from  the  DR  process  missing  the  detector.  A  DR  event  ap¬ 
pears  in  the  pulse-height  spectrum  at  the  total  beam  energy 
unless  one  or  more  of  the  fragments  miss  the  detector.  In  the 
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FIG.  1.  Background-subtracted  pulse-hcighl  spectrum  for  the  DR  reaction  of 
OPCT  recorded  at  =“0  eV  interaction  energy  is  shown  by  solid  squares.  The 
dotted  curve  presents  the  Gaussian  fit  to  the  full  beam  energy  peak. 

latter  case,  the  energy  is  reduced  by  the  percentage  of  total 
mass  which  is  lost.  In  the  pulse-height  spectrum  recorded  for 
the  DR  reaction  of  OPCl+  (Fig  1),  the  presence  of  the  lower 
energy  signals  implies  that  losses  were  occurring  and  the 
contribution  of  these  losses  to  the  total  data  was  calculated  to 
be  16%.  A  Gaussian  best  ht  to  the  total  energy  peak  (Fig.  1) 
would  represent  a  pulse-height  spectrum  in  the  ideal  ease 
when  all  particles  are  detected.  For  cross  section  measure¬ 
ments,  the  loss  of  particles  is  easily  corrected  for  since  it  is 
purely  a  multiplicative  factor.  In  the  case  of  the  heavier 
OPCl2*  ion,  which  has  an  even  lower  kinetic  energy,  the 
measurements  of  both  the  reaction  rate  constant  and  the  re¬ 
action  products  were  even  more  experimentally  complicated, 
since  the  fragment  loss  issue  was  much  more  substantial  than 
for  OPCl+  and  could  not  be  accounted  for  in  the  analysis.  In 
order  to  extract  some  information,  a  detector  with  larger  geo¬ 
metrical  size  was  used.  This  enabled  the  partiele  losses  to  be 
essentially  suppressed,  though,  on  the  other  hand,  this  led  to 
significantly  worse  energy  resolution  due  to  the  increased 
contribution  from  the  electrical  noise  generated  in  the  detec¬ 
tor.  A  Gaussian  fitting  procedure  revealed  that  some  of  the 
fragments  were  still  missed  (Fig.  2),  though  to  a  much  lesser 
extent  than  with  the  smaller  detector.  Moreover,  counts  in  the 
loss  channels  were  dominated  by  the  detector  noise  making  it 
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FIG.  2.  The  solid  squares  show  background-subtracted  pulse-height  spec¬ 
trum  for  the  DR  reaction  of  OPCI2+  recorded  at  *=0  eV  interaction  energy. 
The  dolled  curve  presents  Gaussian  fit  to  the  total  mass  peak  The  abrupt 
cutoff  is  due  to  the  electrical  noise  becoming  comparable  to  the  DR  signals. 

impossible  to  quantify  the  amount  of  loss.  Therefore,  the 
cross  section  measurements  represent  a  lower  limit  to  the 
real  values. 

Branching  fractions  were  obtained  by  inserting  a  grid  in 
front  of  the  detector  with  a  transmission  of 
/= 0.297  ±  0.01 5. 22  Each  particle  traveling  through  the  grid 
has  a  transmission  probability  t  of  making  it  to  the  detector. 
The  DR  events  are  now  recorded  at  the  full  beam  energy 
reduced  by  that  percentage  of  the  total  mass  of  the  fragments 
which  are  stopped  by  the  grid.  For  OPC1*,  the  separation  of 
peaks  in  the  pulse-height  spectrum  was  just  at  resolution 
limit  of  the  detector  and  could  be  used  to  provide  a  unique 
product  determination,  while  in  the  case  of  OPCl2+,  such  a 
situation  did  not  exist  since  the  fragment  peaks  were  below 
the  noise  level  of  the  detector 

In  order  to  address  deficiencies  in  the  literature,  the  ther¬ 
mochemical  properties  of  OPCl2+,  OPCP,  and  PC12  were 
calculated  at  the  G3  level  of  theory  3  using  GAUSSIAN03.  The 
optimized  structures  were  calculated  at  the 
MP2(Full)/6-3 1 G(d)  level  as  part  of  the  G3  compound 
method.  Wave  function  stability  was  verified  at  the 
HF/6-31G(d)  level  for  all  of  the  minimum  energy  structures, 
which  is  the  initial  optimization  step  of  the  G3  method.  The 
G3  0  K  energies  and  298  K  enthalpies  arc  listed  in  Table  I 


TABLE  1.  Energies  al  0  K  and  enthalpies  al  298  K  in  hartrees  calculated  using  G3  theory.  The  ionization 
polentials  have  been  calculated  using  the  G3  energies.  The  standard  heats  of  formation  at  298  K  have  been 
calculated  using  lhe  atomizalion  method.  Literalure  values  are  given  in  parentheses.  See  text  for  details. 


Energy  (0  K) 

Enthalpy  (298  K) 

(kJ  mol"') 

Ionization  potential  (eV) 

opci+ 

-876.080  906 

-876.076452 

827.9  (859s) 

OPCl 

-876.499  870b 

-876.495344b 

-269.5  (— 274,c  -25  ld) 

11.40  (11. 5d) 

OPCl2+ 

1336.221  856 

1336.216403 

550.2 

OPCl2 

-1336.548  600" 

-1336.5429900" 

-303  1 

8.89 

PCL 

-1261.338  130 

-1261.333342 

-76.2 

“Experimental  estimate  (Ref.  25). 
bG3  calculation  (Ref.  8). 


cExperimenial  estimate  (Ref.  24). 
d Experimental  estimate  (Ref  26) 
CG3  calculation  (Ref.  7). 
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FIG  3  Minimum  energy  struclures  for  (a)  0PCI2+  and  (b)  OPCT  optimized 
at  the  MP2(Full)/6*31G(</)  level  as  part  of  G3  energy  calculations. 


for  the  three  molecules  just  mentioned  along  with  previous 
G3  calculations  for  OPCl2  and  OPC1  Refs.  7  and  8  as  well  as 
available  experimental  energetics  for  OPC1  and  OPCF.24  26 
Standard  heats  of  formation  at  298  K  for  the  three  species 
were  calculated  using  the  atomization  method  7  and  are  also 
shown  in  Table  I.  Average  absolute  deviations  of  the  enthal¬ 
pies  of  formation  calculated  using  G3  theory  are 
+3.93  kJ  mol-1  (+0.94  kcal  moP1)23 


III.  RESULTS  AND  DISCUSSION 

The  optimized  geometries  obtained  from  the  G3  calcu 
lations  for  OPCl2+  and  OPCF  are  shown  in  Figs.  3(a)  and 
3(b),  respectively,  having  the  energetics  listed  in  Table  I.  The 
bond  lengths  and  bond  angles  for  neutral  PC12  and  OPCl2 
calculated  by  the  MP2(Full)/6-31G(d)  optimization  step 
agree  within  0.005  A  and  0.5°,  respectively,  with  structures 
calculated  by  Miller  et  al.9  using  the  bigger  6-311+G(d) 
basis  set  and  are  consequently  not  shown  in  Fig.  3.  Similarly, 
the  OPC1  neutral  structure  has  been  calculated  as  part  of 
previous  G3  calculations  8  OPCl7*  is  shown  in  Fig.  3(a)  hav¬ 
ing  a  planar  C2u  structure  that  is  quite  different  from  the 
nonplanar  C5  structure  of  OPCl2.9  However,  the  P-0  double 
bond  in  the  cation  has  essentially  the  same  length  as  in  the 
neutral  molecule,9  while  the  P-Cl  single  bonds  contract  to 
around  1.926  A  in  OPCl2\  In  OPC1+,  the  bond  angle  de¬ 
creases  by  9°  from  the  neutral  OPC1,  while  the  P-Cl  single 
bond  contracts  to  1.948  A  in  the  cation  The  P-0  bond  in 
OPCP  increases  to  an  elongated  double  bond  with  length  of 

I. 566  A. 

The  standard  heats  of  formation  at  298  K  and  adiabatic 
ionization  potentials  (IP)  of  the  neutrals  calculated  using  G3 
theory  are  also  shown  in  Table  I  along  with  the  available 
literature  values.  The  neutral  OPC1  heat  of  formation  found 
using  G3  theory  is  -269.5  kJ  mol-1,  in  agreement  within  the 
uncertainties  of  the  experimental  values  of  -274.15  and 
-251.30  kJmoP1. 24  26  The  calculated  IP  of  OPC1  is 

I I. 40  eV,  in  reasonable  agreement  with  the  1 1.5  eV  value 
estimated  from  electron  impact  experiments.25  This  esti¬ 
mated  value  has  been  used  to  determine  the  OPCF  heat  of 
formation  to  be  859  kJ  mol-1,  which  is  over  30  kJ  moF1 
higher  than  the  G3  value.  However,  the  estimated  value 
agrees  within  the  uncertainties  of  the  experimental  data  used 
to  determine  it.  Previous  comparisons  of  G3  thermochemis¬ 
try  with  the  experimental  values  for  POtClv  species  have 

^  *  7  8  28  29 

typically  shown  decent  agreement.  ’  '  ’  Thus,  the  theoreti¬ 
cal  values  given  in  Table  I  have  been  used  to  determine  the 
energetics  of  the  observed  reaction  channels. 


FIG.  4.  Measured  cross  section  for  the  DR  of  OPCT  is  presented  by  the 
solid  squares  The  line  shows  the  best  fit  [Eq  (4)]. 


The  derived  cross  section  as  a  function  of  the  interaction 
energy  for  the  DR  reaction  of  OPCF  is  shown  in  Fig.  4. 
Since  cross  sections  for  heavy  polyatomic  ions  usually  fol¬ 
low  a  power  law  rr0/E",  the  best  analytical  description  of  the 
data  presented  in  Fig.  4  is  found  to  be 


<t<E)  = 


2.59  X  IQ-16 
E(eV)1  19 


cm2. 


(4) 


The  cross  section  is  described  by  the  same  analytical  expres¬ 
sion  over  a  wide  range  of  interaction  energies,  which  enables 
the  thermal  rate  constant  a  to  be  calculated  analytically  as  a 
function  of  the  electron  temperature  Te, 


a(Te)  =  <r0e" 


T(2  -  n)(kTe) 


1/2-"cm3  s'1. 


(5) 


where  1’  is  the  gamma  function.  Consequently,  taking  into 
account  the  loss  factor  of  16%,  the  reaction  rate  constant  is 
derived  to  be 


/  j  \  -0.89 

a(Tf)  =  7.63  X  10'  (  cm3s~'.  (6) 

As  mentioned  previously,  the  signal  level  for  OPCl2+  was 
close  to  the  noise  level  of  the  detector  (Fig.  2)  which  makes 
it  impossible  to  determine  the  loss  factor.  Assuming  that  all 
of  the  DR  products  are  detected,  the  derived  cross  section  for 
the  DR  reaction  of  OPCl2+  is  shown  in  Fig.  5.  However,  the 
cathode  voltage  needed  to  match  the  velocity  of  the  ion  beam 
was  extremely  low,  4.8  V,  and  the  bit  errors  in  the  DAC 
driving  the  cathode  voltage  give  rise  to  instabilities  in  the 
voltage  and  therefore  uncertainty  in  the  electron  energy.  This 
is  also  the  case  for  OPCF  (the  low  energy  part  of  Fig.  4 
below  3  mcV),  though  to  a  much  lesser  extent  due  to  the 
higher  values  of  the  cathode  voltage.  Consequently,  the  cross 
section  for  the  DR  of  OPCl2+  should  be  considered  as  being 
reliable  only  from  0.006  to  0.3  eV.  The  best  fit  over  this  in¬ 
teraction  energy  range  was  determined  to  be 
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FIG.  5.  The  solid  squares  present  the  cross  section  for  lhe  DR  of  OPCk\ 
The  solid  line  is  Eq.  (7). 


FIG.  6.  Background-sublracled  pulse -height  specirum  for  OPCP  recorded  al 
*=0  eV  interaction  energy  when  the  grid  is  inserted  in  front  of  the  detector  is 
shown  by  solid  squares.  The  dotted  curves  present  Gaussian  curves  of  five 
peaks.  The  solid  curve  is  an  overall  fitting. 


<*£)  = 


6,44  X  IQ-17 
£(eV)' 72 


enr 


(7) 


Simtlarly,  using  Eq.  (5)  with  the  cross  section  presented  in 
Eq.  (7),  the  lower  limit  of  the  rate  constant  for  the  DR  of 
OPCl2+  is  derived  to  be 

a(Te)>  1.2  x  lO-6^)  cm3  s-1 .  (8) 


The  systematic  error  in  the  rate  constants  for  both  OPCF  and 
OPCl2+,  which  is  dominated  by  the  uncertainties  in  the  inter¬ 
action  region  length  and  the  ion  and  electron  currents,  is 
estimated  to  be  =*15%. 

The  DR  reaction  usually  proceeds  through  two  compet¬ 
ing  processes  called  the  direct  and  indirect  mechanisms.  0  31 
The  direct  mechanism  involves  the  formation  of  a  repulsive 
doubly  electronically  excited  state,  while  the  indirect  mecha¬ 
nism  proceeds  through  predissociative  vibrational ly  excited 
Rydberg  states  on  the  neutral  molecule.  It  is  predicted  that 
the  energy  dependence  of  the  cross  section  of  the  purely 
direct  process  is  ME  (Ref.  32)  and,  consequently,  the  tern 
perature  dependence  of  the  rate  constant  must  be  I/7^s  [Eq. 
(5)].  For  the  indirect  mechanism,  this  dependence  can  be  as 
steep  as  l/7^*\31  The  obtained  temperature  dependencies 
[Eqs.  (6)  and  (8)]  imply  that  the  indirect  mechanism  plays  a 
significant  role  in  the  DR  reaction  of  both  OPCl+  and 
OPCl2+.  It  is  noted  that  the  DR  rate  constants  at  300  K  for 
two  other  heavy  triatomic  molecular  ions  03+  and  N3+  are 
reported  to  be  7.37  X  1CT7  and  6.47  X  10“7  cm3  s~!,  respec¬ 
tively,  values  which  are  similar  to  those  obtained  here  for 
OPCl+. 33  34  The  analysis  for  OPCP+  showed  that  the  DR  rate 
constant  is  >  1.2  X  10~6  cm3  s"1  at  an  electron  temperature  of 
300  K.  Such  large  rate  constants  have  previously  been  ob¬ 
tained  for  other  ions  containing  heavy  atoms.  DOCO+,35 
N2OD+,36  DCCCN+37  and  (NO)/.38  The  energy  dependence 
reported  here  for  OPCl2+  is  similar  to  that  found  for 
(NO)2+  3X  The  derived  rate  constants  for  both  OPCP  and 
OPCl2+  indicate  that  the  rate  constant  increases  when  the 
additional  Cl  atom  is  attached  to  OPCP.  Such  a  tendency  has 
previously  been  observed,  in,  e.g.,  oxygen  and  nitrogen  com¬ 


pounds  where  the  rate  constants  increase  in  the  order 
«(X2+)<a(X/)<tt(X4+)  (X=0  or  N).33  34  39  42 

The  background-subtracted  pulse-height  spectrum  for 
OPCP  with  the  grid  inserted  in  front  of  the  detector  is  shown 
in  Fig.  6.  The  spectrum  can  be  separated  into  five  peaks  (Fig. 
6).  Not  every  mass  can  be  identified,  c.g  ,  the  second  peak 
corresponds  to  either  a  single  P  or  Cl  fragment  hitting  the 
detector.  The  branching  fractions  can  be  determined  by  solv¬ 
ing  a  system  of  linear  equations. 


M  O) 

A(P,C1) 

A(0  +  P.O  +  Cl) 

/\(P  +  C1) 

A(P  +  Cl  +  0) 
r  -* 

/(I  -  r)2  /(I-/)  0 

2i(  1  - 1)2  0  r(l  -i) 

=  2i2(  \  -t)  0  /(1-r) 

r2(  1  - 1)  i(\-t)  0 


0 

/( l  -  0 
/(l  - 1 ) 
o 


A7,, 

Nh 

/V, 


,  (9) 


where  Na,  Nh ,  jVc,  and  N(f  represent  the  branching  fractions  of 
the  (la)-(ld)  channels,  A  is  the  number  of  counts  in  the 
peak,  and  /  is  the  transmission  of  the  grid.  It  is  noted  that  the 
third  and  the  fourth  columns  in  the  probability  matrix  are  the 
same,  which  means  that  it  is  not  possible  to  derive  the  frac¬ 
tions  of  the  (lc)  and  (Id)  channels  separately  and,  therefore, 
only  their  sum  can  be  calculated.  In  order  to  determine  the 
branching  fractions  of  all  the  channels  either  the  P  and  Cl 
peaks,  or,  P  +  O  and  Cl+O  peaks  would  need  to  be  resolved, 
something  which  is  not  possible  with  the  current  detector. 
Instead,  we  rely  on  chemical  intuition  based  on  the  limited 
amount  of  data  available  on  DR  studies  involving  triatomic 
systems  with  only  heavy  atoms.  It  is  unlikely  that  during  the 
short  time  period  during  which  recombination  takes  place 
that  the  heavy  outside  atoms  can  combine.  The  three  studied 
systems  confirm  the  absence  of  rearrangement,  i.e.,  02  for¬ 
mation  from  the  DR  of  both  S02+  and  C02+  as  well  as  pro- 
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duction  of  OS  from  the  DR  of  OCS*  have  been  found  to  be 
vanishingly  small. 4^  4<i  Based  on  these  observations,  (Id) 
channel  is  assumed  to  be  negligible.  Finally,  the  branching 
fractions  were  derived  to  be  /Vfl  =  0. 16  ±  0  07,  Nh 
=  0.16 ±0.03,  jVc=0.68 ± 0.05,  and  Nd=0. 

Of  the  two  channels  that  produce  diatomic  neutrals,  the 
more  exothermic  channel  is  more  abundant.  This  result 
might  be  expected  since  the  P-Cl  and  O-P  bond  lengths  are 
1.95  and  1.57  A,  respectively  [Fig.  3(b)].  The  likelihood  of 
the  bond  rupture  in  the  DR  reaction  seems  to  depend  on  the 
bond  strength,  i.e.,  the  cleavage  of  the  weaker  bond  is  more 
likely.  The  branching  fraction  for  the  three-body  channel  is 
16%  which  is  comparable  to  that  of  single  O-P  bond  rupture. 
The  efficiency  at  which  the  DR  reaction  breaks  multiple  co¬ 
valent  bonds  in  light  polyatomic  molecular  ions,  especially 
of  the  type  XH2+  where  X=H,D,0,C,N,P,  etc.,  is  well 
known, 14-1 6,46  and  such  reactivity  has  also  been  observed  in 
heavy  triatomic  ions.  Of  particular  interest  here,  it  has  been 
found  that  the  DR  of  N^+  and  03+  ions  gives  rise  to  almost 
completely  opposite  branching  behavior,  with  94%  of  ozone 
molecules  breaking  both  bonds  in  the  reaction  while  only  8% 
of  azide  molecules  manage  this.3114  Furthermore,  the  full 
fragmentation  channels  in  the  DR  reaction  of  S02+  and  N204 
are  39  and  16%,  respectively.44,47  Due  to  the  sparse  amount 
of  available  information,  so  far  it  is  hard  to  find  similarities 
in  the  breakup  dynamics  of  the  heavier  triatomic  ions. 

IV.  CONCLUSIONS 

The  present  studies  test  the  high  mass  capability  at  CRY¬ 
RING  for  investigating  the  DR  reaction.  For  OPCl+,  with  a 
mass  of  84  amu,  due  to  a  combination  of  the  slowness  of  the 
ions  and  the  poor  phase-space  cooling  =s16%  of  the  DR 
products  missed  the  detector.  Nevertheless,  this  missing  frac¬ 
tion  is  accounted  for  in  the  analysis  and  did  not  present  a 
serious  problem.  However,  for  OPCl2+  (mass  1 17  amu) 
which  has  an  even  lower  kinetic  energy,  such  losses  repre¬ 
sented  a  significantly  greater  problem,  and  in  order  to  make 
some  of  the  measurements  feasible,  a  larger  detector  was 
used.  Although  the  higher  electrical  noise  generated  in  this 
detector  did  not  allow  the  amount  of  losses  to  be  estimated, 
analysis  of  those  particles  which  were  detected  allowed  a 
lower  limit  of  the  rate  constant  to  be  determined.  The  rela¬ 
tively  high  value  of  the  derived  reaction  rate  constant  is  simi¬ 
lar  to  those  obtained  for  other  molecular  ions  of  a  similar 
size  which  could  indicate  that  the  majority  of  the  DR  prod¬ 
ucts  were  detected.  It  is  not  clear  whether  spending  large 
amounts  of  time  focusing  the  ions  could  alleviate  the  beam 
spread  issue  though,  in  practical  terms,  it  is  complicated  by 
the  fact  that  the  electron  cooling  time  is  much  longer  than  the 
lifetime  of  the  ions  in  the  ring.  The  relative  effectiveness  of 
spending  such  a  long  time  cooling  the  beam  also  depends  on 
the  amount  of  ions,  since  there  have  to  be  a  sufficient  num¬ 
ber  remaining  to  make  a  measurement  feasible.  For  the  very 
heavy  ions  discussed  here,  OPCl+  and  OPCl2\  such  ion  cur¬ 
rents  were  not  feasible  and  so  cooling  for  the  required  length 
of  time  would  have  been  impossible.  Additional  constraints 
on  the  reliability  of  the  data  may  also  arise  from  the  low 
kinetic  energy  of  the  electrons  needed  to  match  the  ion  beam 


velocity  for  near  0  eV  relative  collision  energies.  In  deter¬ 
mining  the  product  distributions  from  the  DR  reaction,  these 
measurements  also  push  the  limit  of  what  is  possible.  For 
OPCl2+  it  was  not  possible  to  derive  any  product  information 
at  all.  Even  for  OPCl+,  where  the  smallest  mass  difference 
between  products  is  6  amu,  the  signals  these  masses  gener¬ 
ated  in  the  detector  blended  together,  and  only  a  Gaussian 
analysis  of  the  peaks  combined  with  the  loss  measurements 
and  previously  observed  chemical  precedents  for  the  absence 
of  rearrangement  in  heavy  ions  in  DR  allowed  for  a  complete 
branching  determination  to  be  made. 
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